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Abstract. We present a self-consistent study of coherently coupled two-component Bose-Einstein conden- 
sates. Finite spin-flipping coupling changes the first order demixing phase transition for Bose-Bose mixtures 
to a second order phase transition between an unpolarized and a polarized state. We analise the excitation 
spectrum and the structure factor along the transition for a homogeneous system. We discuss the main 
differences at the transition between a coherent coupled gas and a two-component mixture. We finally 
study the ground state when spin-(in)dependent trapping potentials are added to the system, focusing on 
optical lattices, which give rise to interesting new configurations. 
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1 Introduction 

The physics of multi-component condensates is very rich 
due to the possibility of vector order parameters and the 
presence of different zero-temperature phases. 

In the cold gases context the ability to tune and to 
engineer single- and two-body properties permits many 
implementations of such a system. This toolbox allows 
to address very different and interesting phenomena like 
Andrccv-Bashkin effect [I], persistent currents [2], (in- 
ternal) Joscphson effect 011], Schrodinger-cat- and twin- 
Fock-like states [5JE| , analogues of quantum gravity [7J|51 
[9], spin textures [10], or the most recent and fashionable 
field of non-abclian gauges [TT] , just to cite a few of them. 

In the present work we consider one of the easiest im- 
plementations, namely a 2-component (spinor) condensate 
with an external field that drives the population transfer 
(spin-flipping) between the two atomic levels. In spite of 
the apparent simplicity of the problem, the physics it con- 
tains is very rich, which is reflected in the vast amounts 
of literature generated in the past decades. In this article 
we aim at providing an understanding of the most funda- 
mental features arising from the coupling. Many proper- 
ties of such a system have been already understood and 
addressed in the literature. For the purposes of this arti- 
cle, the most relevant works refer to the structure of the 
ground state and its excitation spectrum [12, 13, 14, 15, 16 . 

Two component condensates are interesting because 
they constitute the generalization of the well known Rabi 
problem of atom optics to interacting extended non-linear 
systems (see, e.g., one of the first experiments |17j ) and 
many properties can be indeed understood using a Bloch 
sphere representation. Under some circumstances it has 
been shown that two component spinors allow for a de- 



scription in terms of a dressed state basis [TMl"M2"fJll2"l"] . 

Moreover the population transfer between the two levels 
turns out to be described by Josephson dynamics, leading 
to what is known as internal Josephson effects (see, e.g., 

m). 

With the increase of atomic species which can be con- 
densed, spinor systems offer new interesting possibilities, 
for instance when species-dependent external potentials 
can be created (see, e.g., the proposals in [2lfl[2"4"n . The 
relation between the latter system and the condensates 
with artificial spin-orbit-like interactions, which is clearly 
built on the physics of two-component spinor condensates, 
would be useful for the field. 

For these reasons in the following we try to give a self- 
consistent and complete description of the system collect- 
ing together a number of results. Some of the results we 
obtain are scattered in literature and proper references 
are provided. New results on the dynamics of these sys- 
tems and on their behaviour in the presence of external 
potentials are presented. 

The structure of the article will be as follows. In Sec. [5] 
we analyse the ground state of the homogeneous system, 
which shows a phase transition between a neutral (GS1) 
and a polarised state (GS2). The static compressibility 
and susceptibility are also studied. The latter has a diver- 
gence at the symmetry breaking point. 

The elementary excitations of the two component sys- 
tem are addressed in Sec. [3] First we write an effective 
quantum hydrodynamic theory to get an insight into the 
structure of modes, Sec. 13.11 The elementary excitations 
are then studied in detail within a Bogoliubov approach 
in Sec. H 

The knowledge of the the spectrum allows us to cal- 
culate in Sec. 2] the density and spin structure factors. 
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Such quantities can be measured experimentally also in 
the trapped case, being related to the local fluctuations of 
the density and of the polarisation. 

The case of trapped gases is discussed in Sec. [5] First 
we address the case where the potentials acting on the 
two species are the same, Sec. 15.11 In Sec. 15.21 we concen- 
trate on the case of an optical lattice acting only on one 
component. 



2 Ground state of homogeneous spinor 
condensates 

We consider a homogeneous spinor condensate whose two 
components a and b interact both via s-wave contact in- 
teractions and via a coherent coupling. Within the mean- 
field framework the system is described by coupled Gross- 
Pitaevskii equations for the spinor components \P a (r,t) 
and #b(r,i) 



at 



at 



h 2 V 2 

2m 

n 2 v 2 

2m 

fi*&„. 



+ V a +g a \V a \ 2 +g ab \V b f 
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(1) 



&b+ 



(2) 



where m is the atomic mass. The contact interaction cou- 
pling constants are given by gi = 4Trh 2 ai/m, with i = 
a, b, ab, where a a and a b are the s-wave scattering lengths 
for components a and b, and a ab that associated to the in- 
teraction between a and b. The term fl introduces a coher- 
ent coupling between the two components, which gives rise 
to phase correlations between the two fluids, in contrast 
to the density-density correlations coming from the in- 
terspecies interaction g ab . Depending on the physical sys- 
tem, this term can have its origin on either a two-photon 
(Raman) process or a direct Rabi coupling between the 
components. For the homogeneous system, the external 
potentials are V a = V b = 0, which is the situation we 
consider in this article except in Sec. [5] 

Due to the flipping term only the total number of par- 
ticles (total density in the uniform system) n = n a + n b is 
conserved. Thus the chemical potential /i is the same for 
both components and the stationary states evolve as 



-ifit 



(3) 



It is convenient to write the spinor components in terms 
of the density n a and the phase <f> a 



J<t>!T 



(4) 



The ground state of the system has been described 
in the literature [T4l[T5l[T8lfT9] . but here we revisit it in- 
troducing a new convenient notation. The ground state is 
given by the values of densities and phases which minimize 



the energy per unit volume 

e(n a , rib) =^g a n 2 a + \^g b n\ + g ab n a n b + 

+ 2\{2\ cos4>y/n a n b - ^(n a + n b ) (5) 

where we have introduced the phase </> = (frba+^o, in terms 
of the phase difference 4> ba — cf> b — <p a and the phase of the 
Rabi coupling, given by fl = \fl\e l ^ n . The configuration 
with minimum energy corresponds to cos <f> = — 1. For fl 
real ((pa = 0, 7r) this means <fi ba = ir for fl > and <p ba = 
for fl < 0; for fl complex, the equilibrium value of (j) ba is 
such that it satisfies <fi ba + <f>n = (2n + l)ir with n E Z. 
Notice that the condition cos</> = +1 can give rise to an 
extremum of the energy [T2lfr3l[T^] , but it will never be the 
global minimum (in fact it is a saddle-point in the energy 
landscape). 

The equilibrium configuration is then characterised by 
the density difference n a — n b . The structure of the ground 
state is better understood in the symmetric case g a = 
g b = g, when the equilibrium solutions must satisfy the 
equation 



.9 - 9ab + 



\fl\ 



(n a - n b ) = 



y/n a n b/ 

This equation admits the following solutions 
(GS1) n a - rib = 0; 



(6) 



(7) 



(GS2) (n a -n b )±=±nJl- 



2\fl\ 



(.9 - 9ab)n 



, (8) 



corresponding to neutral (GS1) and polarised (GS2) ground 
states. Introducing the parameter g ab = g + 2fl/n one 
finds that GS1 (GS2) has the minimum energy provided 
g a b < g a b (g a b > Sab)- In Fig. Q] wc report the polarisa- 
tion of the condensate as a function of the inter-species 
interaction g ab . At the critical point g ab = g ab there exists 
a bifurcation in the ground state solutions, which hints 
at a second-order phase transition. This bifurcation has 
been measured in the experiment reported in Ref. [J and 
has been predicted to be related to the emergence of a 
Schrodinger cat state [5] . It is worth remembering that in 
a 2-state condensate without the flipping term (i.e., n a , 
n b conserved) the condition g ab = g distinguishes a ho- 
mogeneous from a phase separated state (see, e.g., [H]). 
Clearly the presence of the coherent coupling shifts the 
critical value to higher values of interspecies interaction, 
and prevents phase separation by instead creating a po- 
larisation. 

Let us briefly comment here on the relation between 
the scenario above and the internal Josephson effect (see 
also Ref. [H]), which is usually addressed from Eqs. dU 
and @, recognizing fl as the weak coupling. The phase 
we have named GS1 corresponds to a fixed point in the 
Josephson Hamiltonian around which closed orbits ex- 
ist, with vanishing mean polarisation (or population im- 
balance) and a phase difference around 7r, giving rise to 
plasma-like oscillations. Instead, around the fixed point 
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Fig. 1. Different ground states (GSI and GS2) exhibited by 
the two-component spinor system as a function of g a b/§ab- In 

solid: g a = g b ; In dashed: 5g = Q.lg (see text). In dotted: Fi §- 2 - Density (solid) and spin (dashed) static response func- 



unstable solutions. In all cases ft = 0.1 gn. 



tions as a function of the coupling constant of the interspecies 
contact interaction, g ab , above GSI, for O = 0.1 gn. 



GS2 the polarisation is finite, leading to self-trapping dy- 
namics (not a running phase-mode, though, since the mean 
phase difference is locked to 7r). 

The existence of the phase transition can be more 
clearly identified looking at the susceptibility of GSI, which 
within linear response is given by 



X* 



lim ■ 



n b 



V 



g - 9ab + 2\Q\/n ' 



(9) 



and diverges at the phase transition. The compressibility 
on the other hand reads 



,. n a + n b 
Xd = lim 

r)-y0 rj 



9ab 



(10) 



It diverges for g — > —g a b, indicating collapse of the con- 
densate. The quantity rj is the amplitude of a small pertur- 
bation of the energy, characterising the term —r\(n a ± n b ) 
added to Eq. ([5]) , for a density perturbation or a spin per- 
turbation, respectively. The susceptibility and compress- 
ibility are plotted in Fig.[5]as a function of g a b for a system 
in the symmetric GSI. 

Introducing an asymmetry in the interaction poten- 
tial, i.e., g a 7^ gb the ground state always shows a certain 
polarisation and the degeneracy (bifurcation) is removed. 
We take in particular g a = g and gb = g + 5g > g. For 
GSI and g a b <C g a b & small Sg creates a linear polarisation 
in the system given by 



Sn = n a — rib = — 



Sg 



2g- g ab + 2\n\/n + 5g 



(11) 



Obviously for Sg > we always have n a — nb > (see 
dashed line in Fig. [1]). As g a b gab even a small Sg 
strongly polarises the system, and for g a b > g a b the polar- 
isation fastly saturates to the value of GS2 for Sg = 0. On 
the other hand the lower branch becomes a high energy 
minimum and the rib — n a branch becomes a maximum. 



3 Excitation spectrum 

3.1 Hydrodynamic picture: hybridization of spin and 
density modes 

In the present Section we use a perturbative hydrody- 
namic approach in order to get an insight into the nature 
of the excitation spectrum. In the next Section a Bogoli- 
ubov approach is employed for a more quantitative anal- 
ysis. Let us introduce the fluctuation fields II a and <$> a 
for the densities and the phases, respectively. The various 
contributions to the effective hydrodynamic energy func- 
tional for the excitations, Ehd = Eq + V a b, read f33] 



^o=E / 

u J 



Vab 



a—a,b * 



9ab 



Qn 



2m 



2n n 



(12) 



Q 

2ti 



Th 
>>b 



(13) 



where as usual we keep the terms up to second order in the 
fields and we define n = yfn a nb and c a the speeds of sound 
of the two components when the coupling is switched off. 

The analysis is most enlightening for the symmetric 
case g a = g b = g, when (? a = c 2 = gn/m and n = ri/2. 
The energy is obviously diagonalised by the density (or 
in-phasc), £<j = (£ + £i>)/2, and spin-density (or out-of- 
phase), £ s = (£ a — £b)/v2, fluctuation fields, with £ = 
77, </>. The energy functional can be written in terms of 
the new fields as 



Ehd — 



4m. 
h 2 n 
4m 



(V0 d ) 2 

(v</g 2 - 



mcz 



n 2 s + Qnct>l 



(14) 
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where we have introduced the density c ( i and spin c s sound 
speeds (see Eqs. (|T5*|) and (jlBf ). From Eq. (|l"4|) one imme- 
diately sees that while the density sector remains linear 
and gapless, a gap can appear in the spin sector and the 
two branches cross for a certain value of momentum. In- 
deed the equations of motion derived from the previous 
energy functional give the dispersion relations 



with 



C = 



h a g a n a g a bn a b~\ft\ 

-g a n a -K 



9abn a b 



9abn a b ~ \M\ 9abn a b 

-gabflab -gabTlab + \&\ 



-gabnab —gabUab + \&\ 

hb gb"nb 



-gbrib 



(g + g a b)n 

0Jd = c d \k\ = \j — \k\ 



(15) 



u a = c s Vk 2 + 4mO =J( { ~ 9 9ah)n + -\ 
1/ \ 2m m J 

(16) 



where we have defined n a b = ^/n a n b and 
h 2 k 2 

h a = — h 2g a n a + g a bn b - \i, 

Zm 

h 2 k 2 

h = — h 2g b n b + g ab n a - fJ,, 



2m 



allowing for a value ko at which uj a ' = w s . From the dis- 
persion relations Eqs. (fT5")) and (JTHJ) it is clear that such 
an unpolarised phase can be stable against collapse only if 
g a b > —g and against polarisation only if g a b < g + 2fi/n. 
The latter is precisely the condition we found in the previ- 
ous Section for GS1 to be the lowest energy state. Notice 
that since the leading term at small momentum for the 
spin spectrum is k 2 its coefficient is wrong within hydro- 
dynamics. To get it right it would be enough to include 
the next order term, i.e. (V77 s ) 2 / (2nm) in the energy 
functional Eq. (|T3j), 

Any small difference between g a and g b would couple 
the density and spin-density leading to an hybridization 
of the modes and thus to an avoided crossing between the 
new branches. 



3.2 Bogoliubov excitations 

In this Section we derive the dispersion law of the two- 
component spinor BEC for the excitations above the dif- 
ferent ground states within Bogoliubov approach. The dy- 
namics of the system is given by the coupled time-dependent 
Gross-Pitaevskii Eqs. (|TJ) and ([2]). To address small am- 
plitude excitations above the ground state we write the 
time-dependent wave functions as 

9„ (r, t) = e-^' h [% + e**'6^[T, t)] (17) 

with /i the chemical potential of the complete system, i\) a 
the ground state wave functions defined in Eq. and 
a = a,b. Since we are dealing with a homogeneous system 
we can write the excitations SW^ in terms of plane waves 
with amplitudes u a and v a , 



M= 2 



g a n a + gbm + g a bn - 



\/n a n b 



-h b 
(20) 



(21) 
(22) 
(23) 



We moreover choose the normalisation such as that the 
amplitudes u a and v a satisfy 



KI 2 -kl 2 + kl 2 -M 2 = ±i 



(24) 



3.2.1 Excitations above the symmetric GS1 

Let us first analyse the excitations above GS1 for g a = 
g b = g, which corresponds to n a = n b = n/2 (see Fig. [T|). 
This solution only exists for g ab < g + 2\Q\/n. The disper- 
sion relations are 



(^i) 2 = 

(hu 2 ) 2 - 



h 2 k 2 (h 2 k 



2m 
h 2 k 2 



2m 

h 2 k 2 



+ (g + gab)n 
+ (g - g ab )n + 4\n\ 



2m \ 2m 
+ 2\n\[(g-g ab )n + 2\n\] 



(25) 



(26) 



These results arc in agreement with Refs. p~2lfT3iri4] and 
are plotted in Fig. [3] The frequency u>i corresponds to 
a density mode as can be seen in the bottom-left panel 
of Fig. |3] (see also text below). At low k the dispersion 
is linear with sound speed Cd- Notice that the frequency 
LO\ does not depend on the coupling Q and is given by 
the same expression as the density mode of the uncoupled 
two-component case [5S] (see Eq. (|28|l below). In contrast, 
frequency u>2 - which corrects the expression Eq. (|16[) - 
corresponds to a spin mode (see bottom-right panel of Fig. 
[3]and text below), that goes as ~ k 2 for k — > and has a 
gap at k = of value 



6V„ 



JQsl-t— uit) 



89* = UfT e l ( kr -" t ). 



(18) 
(19) 



frwj = J2\Q\[(g-g ab )n + 2\n\] , 



(27) 



This gap can be shown to correspond to the Josephson 
The equations for the Bogoliubov coefficients u's and v's frequency for small amplitude oscillations. In the limit of 



can be written as 



f Ha ) 




/ Ua 




-< 


v a 


U b 




Ub 


\VbJ 







\Q\ — > we recover the result for a mixture of two com- 
ponents, 



2J.2 rfc2;„2 



(M 2 



K 2 k 2 
2m 



K 2 k 2 
2m 



n(g ± g a b) 



(28) 
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Fig. 3. Top: Dispersion relation e(k) above GS1, for g a b = 
0.5 g and .!? = 0.1 gn. Bottom: Components of the eigenvectors 
corresponding to eigenvalues u)\ (left) and U2 (right). 
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Fig. 4. Top: Dispersion relation e(k) at the critical point, for 
g a b = g and O — 0.1 gn. Bottom: Components of the eigen- 
vectors corresponding to eigenvalues toi (left) and u>2 (right). 



As already anticipated in Sec. 13.11 there is a crossing be- 
tween the two frequencies at a finite value of fc, 

Notice that kg exists provided g a h < g + 2\f2\/n, or g a b — 
2\Q\/n < g, which is cxaclty the condition for GS1 to be 
the ground state. As the critical condition is approached 
the crossing occurs at lower k and the gap energy ap- 
proaches zero. At the critical condition the gap wj closes 
and the dispersion relation becomes linear at low fc, as can 
be seen in Fig.0] Such a behavior is very different from the 
softening of the mode in the mixture case at the demixing 
point, for which the frequency goes as w ~ k 2 at low k. 
Finally, for g a b > g a b the frequency cu 2 (calculated above 
GS1) becomes imaginary, leading to instability, since the 
real ground state under this condition is GS2. 

To understand the density and spin characters of the 
modes u>\ and ui 2 it is useful to look at their eigenvectors 
(bottom panels of Fig. [3]). The eigenvector corresponding 
to <jj\ is a density-like (in-phase) excitation which satisfies 
u b = u a = Uk and Vb = v a = Vk and shows the usual 
infrared divergence — > fc -1 / 2 , Uk —> —v^ for k 0. The 
eigenvector corresponding to ui2 is spin-like (out-of-phase) 
and satisfies Ub = —u a = — Uk and Vb = —v a = —u^. In 
this case, away from the phase transition, the eigenvector 
components are not strongly dependent on k and they 
remain finite over all /c-space at the value given by the gap. 
Approaching the phase transition the spin mode diverges 



for k — > (see bottom left panel of Fig. 0]), but it's and v's 
with the same index a, b have the same sign, which will 
lead to a diverging spin static structure factor (see Sec.|3]). 
It is interesting to note that there is a change of sign in 
u a and Ub for the spin mode as soon as g a b > 9 (compare, 
for instance, Figs. [5] and HJ). 

3.2.2 Excitations above the asymmetric GS1 and GS2 

From the point of view of the excitation modes, the GS2 
phase and the case with g a ^ gb are qualitatively equiva- 
lent since both situations are characterised by a finite po- 
larisation and a coupling between density and spin-density 
modes (hybridization). Full analytical expressions for the 
frequencies oj\ and lo 2 in the general case were reported in 
Rcf. [14] . but here we calculate them numerically to have 
a clearer access to the physics. An example of the result- 
ing avoided crossing of the modes uj\ and 101 is given in 
Fig. [5] for the case gb — g a - The fact that 8g > shifts 
the value of k at which the avoided crossing occurs at 
k > ko- For 5g < 0, the avoided crossing takes place at 
k < kg. The bottom panels of the figure show the com- 
ponents of the eigenvectors corresponding to the new lu\ 
and UJ2 modes. We can see that away from fco their val- 
ues approach (perturbatively) the values of the symmetric 
case (Fig. , while around fco they show a sharp change 
of magnitude. This change is steeper as 8g — > 0, becoming 
a delta function for 5g = 0. Notice that this fact is hidden 
in Fig. [3] into the definition of the modes u>\ and u>2 before 
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Fig. 6. Dispersion relation e(k) above GS2 for g ab = lAg and 
J? = 0.1 gn. The polarisation for this case is (n a — nt)/n — 
—0.866. Bottom: Components of the eigenvectors correspond- 
ing to eigenvalues wi (left) and 012 (right). 



the symmetric GS1 case vl is given by 



and after the crossing (compare solid and dashed lines in 
Figs. Eland EJ. 

For the GS2 the point where the avoided crossing oc- 
curs depends strongly on the specific location along the 
bifurcation curve. For n a — small (that is for g a b > g a b) 
the avoided crossing will show up at low k and the energy 
difference between the two modes will be small. Instead, 
for a polarisation tending to unity, the avoided crossing 
moves to k — s- oo and the energy difference between the 
two modes increases. We can see from the bottom panels 
of Fig. |5] that when the polarisation is large (and the en- 
ergy gap between the two modes at the avoided crossing 
grows) the behaviour of the eigenmodes departs from their 
behaviour in GS1. 



3.2.3 Landau critical velocity 



v L =< —y/2\{2\[(g- 9a}) )n + 2\{2\] + 



1/2 



— [{g-9ab)n + A\Q\ 
2m 



(30) 



In the limit SI — s- we recover the sound velocity of the 
spin mode for a mixture of two condensates, see Eq. ((281) . 
For g ab g ab we get v L -> (g ab - g)n/m. Notice that 
while Vl =/= at the transition point at finite f2, for a 
mixture (or SI = 0) vl =0. 

Clearly a pure density probe would not be very effec- 
tive to measure the reduction of vl since it comes mainly 
from spin-like excitations. One could think to use a probe 
that acts only on one of the two components in such a 
way to be well coupled to both the density and the spin 
sector. Recently this topic has received some attention in 
a two-component system with spin-coupling in a z [27] . 



One of the consequences of the spectrum discussed above 
is that the Landau critical velocity vl — min[cj(fc)/fc] 
is not just the slope of the density mode, namely vl ^ 
linifc-^o tJi(k)/k = Cd, but it is smaller. In particular for 



4 Structure factors 

The dynamic response of a system to an external per- 
turbation can be studied within linear response theory in 
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terms of the dynamic structure factor 

S F (k,u) = \Y. \{™\F{k)\0)\ 2 5{u-Lo nM ) (31) 

where F(k) is the operator related to the perturbation one 
is interested in (see for instance Refs. [28ll29j ) . Integration 
of Sp(k,ijj) over the frequency w gives the static structure 
factor Sp{k), which quantifies how strongly the different 
modes are excited by the perturbation. 

In a two-component spinor condensate the interesting 
response functions are related to the in-phase and out-of- 
phase density modes, whose corresponding operators are 



F d (k)=5p\k) = 



E 

(T—a.b 



(32) 



F s (k) = a z 5p f (k) = ^2 sg 11 ^) (aL«o<r + aL^-far 



(T—a,b 



(33) 



a, 



where a ka (a^) creates (annihilates) a particle a 
in a state with momentum k, a z — {{1 0}, {0 — 1}} is the 
third Pauli matrix, and sgn(a) = — sgn(6) = 1. Within the 
Bogoliubov approximation discussed in Sec. 13.21 we can 
introduce the quasi-particlc creation (annihilation) oper- 



ators, bl (bka) for the modes a = 1, 2 



"ka 



E 

a=l,2 



-i<t><7 



a=l,2 



{Ukaaha + V*_ k(7a b f _ ka ) 
V-kaab-ka) 



£t 

kaa u ka 



(34) 
(35) 



where {iikaon Vkaa} are the eigenvectors related to the 
a mode found in Sec. 13.21 The density and spin dynamic 
structure factors are easily calculated and can be written, 
respectively, as 



S d (k,u) =S d (k)S(uj - Wl ) + S 2 d (k)S(uj - oj 2 ) 
S s (k,cj) =Sl(k)S(uj - wi) + S 2 (k)S(uj - uj 2 ) 



with 



S2(k) = 



Sf(k) 



^ (Ukaa + Vkaa ) 

(7 



(36) 
(37) 



(38) 



(39) 



We immediately see that for the symmetric GS1, due to 
the symmetry of Bogoliubov amplitudes u and v, the den- 
sity structure factor is completely exhausted by the den- 
sity mode u)\ , while the spin structure factor is completely 
exhausted by the spin mode ui2- On the contrary both for 
GS2 and the asymmetric case (i.e., g a ^ gb), the inter- 
ference term in the equations becomes important due to 
mode hybridization. 

The static structure factors are simply 



S d ( s )(k) 



(40) 




1.2 
1.0 
0.8 
0.6 
0.4 
0.2 
0, 





(c) 






- 














Fig. 7. Density (continuous black line) and spin (red dashed 
line) static structure factors for g a — gb for GS1 with (a) 
gab/g a = 0.5, (b) g a b/g a = 1.15 and for GS2 with (c) 
gab/ga = 1.25, (d) g a b/g a = 1.5. In all cases, Q = 0.1g a n 
and gab/ga = 1.2. 

and are plotted for g a = g\> in Fig. [7]as a function of k, for 
four different values of g a b before and after the transition. 
As expected the density structure factor shows the usual 
linear behaviour at low k found for a condensate, and it 
satisfies Feynman criterion (see for instance Ref. [2"5]). 

The spin structure factor, instead, shows a gap at k = 
which arises from single-particle coherence (see also next 
section). This behaviour of the spin structure factor is 
fundamentally different to what is usually found in Bosc 
gases at T = 0, and in particular in mixtures of two com- 
ponents, where the spin structure factor fulfills Feynman 
rule (see for instance 30 ). At g = g a b the structure fac- 
tor reaches unity and keeps constant for all k. This is 
related to the change of sign in u a and Ub for g = g a b (see 
Sec. I3.2.T]) . Approaching the phase transition, Fig. jT^b), 
the spin structure factor presents a maximum for k — > 0, 
which diverges at precisely g a b = (jab- This is due to the 
fact that close to the phase transition spin fluctuations 
are very large and, since they are proportional to the spin 
static structure factor (see next section), this latter is also 
very large. Above the transition point, Fig.fTJc), the spin 
structure factor maximum moves toward k ^ due to the 
polarisation of GS2. Finally when the system lies deep 
in the GS2, Fig. [7] (d), the spin-mode reduces essentially 
to a single-component density mode and for this reason 
behaves linearly at small k and the gap tends to zero. 

The effect of finite temperature T well below the crit- 
ical condensation temperature T c can be introduced by 
using detailed balanced. The static structure factors read 



Huj a (k) 



S v {k,T)= J2 Smooth- 

a=l,2 B 



d,s. (41) 



The zero-temperature results are recovered for ksT <C 
huj a (k). For large Q (deep in GS1) the latter condition 
can be easily satisfied for the spin-mode. This is not true 
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close to phase transition. On the other side of the phase 
transition the maximum at finite k will be hidden by the 
usual enhancement of the static structure factor for k — > 0. 



4.1 Density and polarisation fluctuations 

The peculiar behavior of the structure factors will show 
up in the local density AN 2 and polarisation AM 2 fluc- 
tuations of the spinor gas. Such fluctuations, which are 
nowadays experimentally accessible (see e.g., [511152,33, 
34,35,36 ), arc important quantities characterising a sys- 
tem (for a recent review in the context of cold gases see 
[37]). If R is the linear size of the spot where the fluctua- 
tions are measured one can approximate the fluctuations 
as 

AN 2 = nJ S d (k,T)H(k)j— ~ NS d (l/R,T),(42) 

/elk 
S s (k,T)#(k)— ^ ~ NS S {1/R,T);(A3) 

with D the dimensionality of the system and H(k) the 
Fourier transform of a geometrical factor that depends 
on the shape of the probe cell where the measurement 
is performed. The last approximate equality [57], where 
R is the linear size of the spot, has been recently used 
to probe the (density) structure factor of a 2-dimensional 
Bose gas [35] . For a spinor gas the local fluctations of the 
polarisation can be also used to measure the spin structure 
factor. 

The behaviour of S s (0) can be easily understood by 
looking at the spin fluctuations. Indeed, if we look at them 
at the single-particle level (two level system), we always 
find that AM 2 = (a 2 ) - (a z ) 2 ^ 0. Since we have the 
relation S s (0) = AM 2 it follows that the spin structure 
factor is gapped. In this sense we can say that the gap in 
S s comes from single-particle coherence. 

On the other hand since the phase transition is due 
to spin-density instabilities, as we have already seen from 
the susceptibility, Sec. [2] close to the phase transition the 
fluctuations of the polarisation, i.e., S s (0), diverge. For 
instance, in the context of Bose-Bose mixtures such fluc- 
tuation enhancement has been very recently observed in 
a quenching experiment |38j . 



5 Trapped two-component spinor condensates 



In this section we discuss how the GS1-GS2 picture pre- 
sented in Sec. [2] can be applied to inhomogeneous sys- 
tems. We find that local density approximation captures 
the main features of the system brought about by the 
trapping potentials. Within this approach the equations 
describing the ground state of the trapped system can be 
found from Eqs. JT]) and @. Neglecting the kinetic energy 



we find by addition and subtraction the two equations 



\o\ 



g - g a b 



(n a - Tib) =Vb~ V a , 



(44) 



s]n a n hi 

g + g ab - -pd= ) K + ri b ) = 2p - (V b + V a ). (45) 
\/n a n b 



For V a — Vb — we obviously recover the results from 
Sec. [U To exemplify what happens for V a ,Vb ^ 0, we 
consider two situations. In the first one, the trapping is 
harmonic and equal for species a and b. In the second 
example, instead, one of the components is subject to an 
optical lattice potential, which is not directly seen by the 
second component. In both cases, numerical solutions of 
the full GP equations are provided for comparison. 



5.1 Harmonically trapped two-component spinor 

Let us consider what happens if there is a spherically sym- 
metric external harmonic confinement Vb (r) = moj 2 r /2, 
which acts in the same way on both spinor components. 
We first address this problem within local density approx- 
imation. For V a = Vb the right-hand side of Eq. (|4"4"]) van- 
ishes, which means that there are two possible ground 
states characterised locally by the same solutions as the 
homogeneous system, Eqs. ([7J and ©. This allows us to 
introduce the local critical value g ab (r) = g + 2\Q\/n(r), 
with n(r) the total local density. Analogously to the ho- 
mogeneous case, for g ab {f) > g a b the system is locally in 
GS1, while for the opposite condition the system in locally 
in GS2. Since at low density the Rabi flopping term always 
dominates, when a two-component spinor condensate is in 
a trap two possible scenarios can exist: 

(i) the whole system is in GS1, i.e., the critical condi- 
tion g a b < g ab (0) is fulfilled; 

(ii) GS2 is the lowest energy state in the center of the 
trap, i.e., g a b > g ab (0). Then there always exsists some 
critical radius R c fulfilling that the (decreasing) density is 
such that g ab < gab{r) for r > R c and GS1 is the lowest 
energy state. In this case there is coexistence of the two 
phases and the critical radius is given by g ab = g a b{Rc)- 

The first scenario is similar to the usual Thomas-Fermi 
approximation for a single condensate. We concentrate 
therefore on the second situation. For large r the system 
will be in GS1, that is, n a (r) = rib(r) = n(r)/2. For a 
spherically symmetric harmonic potential the density in 
the GS1 phase can be calculated from Eq. (|4"5|) as 



n(r) 



M - Vho(r) 



g + g ab - 2\n\/n(r) 



(46) 



At the critical radius the density fulfills n(R c ) = 2| Q\/(g a b- 
(?) and substituting above we find 



R c 



2^ + 4|/2| 
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g a b 



g > g a b- 



(47) 



For r > R c the system is in GS1, while for r < R c GS2 
is the ground state and the system is polarised. Figure [5] 
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Fig. 8. Example of the ground state in a harmonic trap, with 
an inner region described by GS2 and an outer region described 
by GS1. Parameters: g a t/g — 1-3, (2/fi = 0.1 and R c /Rtf — 
0.55. The quantities no and Rtf represent, respectively, the 
density at r = and the radius where the density vanishes. 




Fig. 9. Two-phase ground state density profiles of compo- 
nents a and b from the numerical solution of one-dimensional 
GP equation. The parameters are: g a t/g = 1.25, fi/fi = 0.05 
and R c /Rtf = 0.55. The quantities no and Rtf represent, 
respectively, the density at x = and the position at which 
the density vanishes. 



shows the density distribution of components a and b, as 
well as the total density (in the inset), of a 2-phase config- 
uration. Due to the use of the local density approach the 
density profiles measurement is a direct mapping of the 
phase diagram of the homogeneous system, Fig. [TJ 

Including quantum pressure the density profiles are 
smoothed out. This behavior is shown in Fig. El obtained 
by numerically solving the quasi-onc-dimcnsional spinorial 
GP equations in a trap. We have checked that equivalent 
density distributions are found in the 3D case for a spher- 
ically symmetric trap. 

The two-phase configuration that we find is reminis- 
cent of the results reported in Ref. |39j in the intermedi- 



ate state in the transition from the misciblc to immiscible 
regimes (see also the experiment of Ref. [21]). 



5.2 Spinor condensate in a spin-selective optical lattice 

Interesting scenarios appear when the two components feel 
different external potentials. We concentrate on the situ- 
ation where component a feels an optical lattice poten- 
tial, while component b does not. Since the atomic transi- 
tion is locally out-of-resonance, polarised ground states 
are favoured, as can be seen from Eq. (|44[) . Assuming 
| SV | = 114 — 141 <C |14 + Vi,| we can perform a perturbative 
analysis. We define V a = Vq + 8V C , with a = a, b, which 
allows us to write the densities as n a = n£ + 8n a , where 
the n° are the unperturbed densities, that is the densities 
corresponding to Vq. For GS1 n°(x) = n°(x) = n (x)/2 
and we find at first order 



8n<j 



g/x - AV/2 

9 + 9ab 



sgn(cr) 



SV/2 



2\S7\/n Q + g - g ab ' 



(48) 



where Sfi is the change in chemical potential and AV = 
8V a + SVb. Notice that from this expression the external 
potential acts on component a but also on component b. 
This means that the interspecies interaction as well as the 
Rabi coupling create an effective external potential for the 
untrapped component. 

The first term in the right-hand side of Eq. (|48[) acts as 
renormalized chemical potential that takes into account a 
mean potential change, AV/2, and gives the same contri- 
bution to the densities. The second term instead is species- 
dependent, showing a positive or negative contribution 
depending on the component. When the first term domi- 
nates, the fluctuations 8n a have the same sign, while when 
the second term dominates 8n a and 5rib are out-of-phase. 
It is clear then that for large fl or cquivalently small den- 
sities the system tends to be unpolarised; in contrast, for 
high densities (for instance, at the trap center) the system 
tends to be polarised. 

Obviously the density change in component a is oppo- 
site to 8V and close to the phase transition even a small 
difference in the external potential gives rise to a large lo- 
cal polarisation, and eventually the perturbative approach 
breaks down. A similar solution and argument hold for 
GS2. 

To test the predictions of local density approximation, 
we numerically find the stationary ground state of the 
(quasi-one-dimensional) spinor Gross-Pitaevskii Eqs. (TT])- 
J2]) for potentials V a = Vho + V Q i and Vb 
V i(x) = V (sin 2 (k L x) - 1/2) and V ho (x) 
In this case it is clear that for small g a b the system is 
everywhere in the state GS1 and the effect of the Rabi 
coupling is just to make the state b to feel the lattice, 
obtaining a polarised in-phase periodic spinor condensate, 
as the one shown in Fig. [TO] On the other hand if g a b is 
large enough we have a central region where the system is 
in the GS2 state, which results in a polarised out-of-phase 
spinor condensate as shown in Fig. [11] As discussed above 
when the density is small the system enters in GS1 and 



Vho, with 
moj 2 x 2 /2. 
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Fig. 10. Density profiles (solid line) for a spinor condensate in 
a spin-dependent optical lattice plus harmonic trap (see text) 
in GS1. The parameters are g a b = 0, Q/(gn) = 1, Vo/(gn) = 1 
and d = 0.8 a,ho, with a; 10 = ^h/mw the oscillator length. 
Dark line shows n a and light (red) line shows rib- We report 
for comparison also the case with Q = (dashed lines). 




Fig. 11. The same as in Fig. 1101 but in coexistence case of both 
GS1 and GS2. The parameters are : g a b/g = 40, Vo/(gn) = 1, 
Q/(gri) = 1 and d = 1.6 dho, with a; 10 = yjh/mui the oscillator 
length. Dark (black) solid line shows n a and light (red) solid 
line shows rib- 



this is clearly seen in the tails of the condensate, where 
the oscillations of the two components are in-phase. 



complete description and interpretation of the fundamen- 
tal properties of the system. 

In light of our results, we can compare a 2-componcnt 
condensate with and without the Rabi coupling term Qa x . 
In the former case the transition is second order leading 
to a partial global polarization, while in the latter case it 
is first order and coexistence of domains of polarization 
±1 are present in the system. The spectrum shows a gap 
for Q ^ since it breaks the 0(2) symmetry transverse to 
the polarization, making excitation about a x (a y ) massive 
(Fig. [3]). The second order phase transition is due to the 
closing of the spin gap when the inter-species interaction 
is large enough (see Fig. 0]). For Q = the spectrum is 
gapless since any angle in the a x -a v plane (of the Bloch 
sphere) is equivalent. In this case the phase transition is 
driven by the softening of the phonon-likc spin-mode. 

The excitation spectrum behaviour is reflected in the 
density and spin static structure factors. The first one has 
the same behaviour as in a single component condensate 
and satisfies the Feynman relation. The spin structure fac- 
tor instead shows a peculiar behaviour and in general does 
not follow the Feynman relation. In particular due to the 
coherent coupling it remains finite at zero momentum, and 
it shows a divergence approaching the phase transition 
(Fig. [7]). Such a divergence is easily explained in terms of 
the fluctuations of the polarization (see Sec. 14.11) . 

We also discuss the coexistence of the two ground states 
once the gas is confined in a harmonic potential. In this 
case the phase diagram can be almost directly read out 
by measuring the local density of the two components, as 
clearly seen in Fig. [8] 

Even more interesting is the effect of a spin-dependent 
external potential, which creates an entanglement between 
the internal and the external degrees of freedom of the gas. 
The possible configurations are many and we concentrate 
on the case of an optical lattice potential Vol felt by one 
(component a) of the two components, i.e., a potential 
Vol <& (1 + cr z ). The results obtained numerically solving 
coupled GP equations are reported in Figs. ITTJlfTTl where 
we show how depending on the regime the periodic struc- 
ture induced on the "free" component b can be either in- 
or out-of-phase with respect to the trapped component a. 

We acknowledge very useful discussions with I. Caru- 
sotto, Y. Li, D. Papoular, S. Stringari and A. Trombcttoni. 
This work has been financially supported by ERC through 
the QGBE grant and by Provincia Autonoma di Trento. 
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